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ABSTRACT 

A  multi-platform  angle-only  tracking  system  combines  the  angular  mea¬ 
surements  from  distributed  networked  sensors  in  order  to  estimate  the  full 
kinematic  target  state.  This  report  investigates  the  effect  of  communication 
bandwidth  on  track  quality  in  a  system  which  receives  angle-only  measure¬ 
ments  from  two  networked  sensors  installed  on  two  airborne  moving  platforms. 
The  investigation  is  based  on  the  Cramer- Rao  lower  bound  of  the  mean-square 
range  error  for  the  case  considered.  The  bound  is  derived  for  recursive  estima¬ 
tors  with  prior  information  and  compared  with  two  algorithms:  (i)  maximum 
likelihood  estimation  over  a  cumulative  measurement  set  and  (ii)  extended 
Kalman  filter  with  triangulated  range  estimates. 
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The  Influence  of  Communication  Bandwidth  on  Target 
Tracking  with  Angle  Only  Measurements  from  Two 

Platforms 


EXECUTIVE  SUMMARY 

Tracking  targets  using  line-of-sight  angular  measurements  is  very  important  for  tactical 
air  surveillance.  One  reason  is  that  passive  sensors  (such  as  the  ESM  receiver,  Infrared 
Search  and  Track  sensors  or  the  radar  in  a  passive  mode),  which  enable  covert  mode 
operation,  provide  only  the  target  angular  measurements.  However,  even  the  use  of  an 
active  sensor  such  as  the  radar,  in  the  presence  of  Electronic  Counter  Measures  (noise 
or  deceptive  jamming)  provides  reliably  only  the  measurements  of  the  target  line-of-sight 
angle. 

The  aim  of  a  multi-platform  angle-only  tracking  system  is  to  combine  the  angular 
measurements  from  (possibly  moving)  distributed  networked  sensors  in  order  to  estimate 
the  full  kinematic  target  state  (including  range,  speed  and  heading).  The  quality  of  target 
tracking  is  obviously  going  to  be  influenced  by  the  amount  of  information  that  is  exchanged 
between  the  networked  sensors.  In  practical  systems  this  amount  of  information  is  limited 
by  the  capacity  of  the  communication  link  or  its  bandwidth. 

This  report  investigates  the  effect  of  bandwidth  on  the  track  quality  in  a  system  which 
receives  angle-only  measurements  from  two  networked  sensors  installed  on  two  airborne 
moving  platforms.  The  method  of  investigation  is  based  on  the  Cramer-Rao  lower  bounds 
which  are  derived  for  recursive  estimators  with  prior  information.  As  expected  the  track¬ 
ing  performance  in  the  limited  bandwidth  environment  is  between  the  two  extremes,  the 
full-connection  and  the  no-connection  cases,  and  it  follows  the  pattern  of  the  received 
measurements  over  the  link.  As  time  progresses,  the  limited  (but  non-zero)  bandwidth 
cases  tend  to  approach  the  full  connection  performance.  This  effect  is  more  pronounced 
when  angular  accuracy  is  higher. 

The  analysis  shows  that  only  relatively  modest  communication  bandwidth,  well  within 
the  capacity  of  modern  data-links  such  as  Link-16,  is  required  for  the  target  state  estimate 
to  rapidly  approach  the  limit  posed  by  the  Cramer-Rao  bound.  Thus  it  appears  that 
netting  of  angle-only  information  could  provide  considerable  tactical  advantage  in  resolving 
target  position,  heading  and  speed  even  in  the  presence  of  ECMs.  Further  investigation 
of  this  research  area  is  warranted  with  an  aim  to  relax  some  of  the  assumptions  made  in 
this  analysis. 
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1  Introduction 


Modern  tactical  surveillance  systems  increasingly  tend  to  operate  in  a  covert  mode 
using  passive  instead  of  active  sensors.  For  air  surveillance  typical  passive  sensors  include 
ESM  receivers,  Infrared  Search  and  Track  sensors  and  the  radar  in  a  passive  mode.  The 
main  disadvantage  of  passive  sensors,  however,  is  that  they  are  unable  to  measure  target 
range;  they  provide  only  target  angular  measurements. 

The  aim  of  multi-platform  angle-only  tracking  system  is  to  combine  the  angular  mea¬ 
surements  from  (possibly  moving)  distributed  networked  sensors  in  order  to  estimate  the 
full  kinematic  target  state  (from  which  one  can  then  work  out  the  range,  heading,  speed). 
The  single  platform  case  has  been  studied  in  numerous  publications,  some  more  recent 
ones  being  [1],  [3,  Ch.5],  [4].  It  is  well  known  that  in  the  single  observer  case  the  ownship 
must  outmanoeuvre”  the  target  in  order  to  estimate  the  full  target  kinematic  state  [5].  In 
the  multi-platform  case  the  ownship  manoeuvre  is  not  necessary  (providing  the  vergence 
angle  is  large  enough)  but  the  quality  of  target  tracking  is  influenced  by  the  amount  of 
information  that  is  exchanged  between  the  networked  sensors.  In  practical  systems  this 
amount  of  information  is  limited  by  the  available  capacity  (or  bandwidth)  of  the  tactical 
data  link  used  for  communication. 

This  report  investigates  the  effect  of  bandwidth  on  the  track  quality  in  a  system  which 
receives  angle-only  measurements  from  two  networked  sensors  installed  on  two  airborne 
moving  platforms.  The  investigation  is  based  on  the  study  of  the  Cramer-Rao  lower  bound 
(CRLB)  which  defines  the  best  achievable  performance.  The  CRLB  in  the  case  considered 
is  shown  to  depend  on  the  observer-target  geometry,  the  angular  measurement  accuracy, 
the  sampling  interval  and  the  bandwidth.  The  report  follows  the  style  and  extends  the 
results  recently  reported  in  [4],  The  bound  is  derived  for  recursive  estimators  with  prior 
information  and  compared  with  two  algorithms:  (i)  maximum  likelihood  estimation  over 
a  cumulative  measurement  set  and  (ii)  extended  Kalman  filter  with  triangulated  range 
estimates.  In  order  to  make  the  analysis  tractable  the  following  assumptions  are  made 
in  the  report:  measurements  are  collected  synchronously  on  both  platforms  with  a  unity 
probability  of  detection  and  with  no  false  alarms;  only  a  single  non-manoeuvering  target 
is  being  tracked;  there  are  no  registration  errors  and  there  is  no  transmission  delay. 

The  report  is  organised  as  follows.  Section  2  defines  the  problem  of  angle-only  tracking 
with  two  communicating  airborne  platforms.  Section  3  derives  the  Cramer-Rao  lower 
bounds  for  the  described  problems.  The  influence  of  bandwidth  on  tracking  performance 
is  studied  in  Section  4,  using  the  derived  Cramer-Rao  bounds.  This  section  also  compares 
two  tracking  algorithms  against  the  bounds.  Section  5  summarises  the  main  results  of  the 
report. 


2  Problem  Description 

Target  T  in  Figure  1  is  moving  with  a  constant  velocity  along  a  straight  line  in  x- 
y  plane,  while  two  friendly  platforms  Pi  and  P2,  connected  by  a  tactical  data  link,  are 
collecting  angular  (azimuthal)  measurements  z\  (k)  and  z2  (k)  respectively  at  time  instances 
4  =  kT0,  k  =  1 , . . . ,  IV.  Here  Ta  is  the  sampling  interval  while  N  is  the  number  of 
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measurements.  Azimuth  measurements  are  defined  with  respect  to  the  y  axis  direction 
and  are  affected  by  additive  noise,  i.e. 

Zm{k)  =  arctan  +  wm(k)  (m  =  1,  2)  (1) 

where  {xt{k)fyt(k)}  are  the  unknown  target  coordinates  and  {xm(k),  ym(k)}  are  the  known 
coordinates  of  sensors  m  —  1,2.  Noise  processes  w\(k)  and  w2 {k)  are  mutually  uncorre¬ 
lated  zero-mean  white  Gaussian  with  variances  a\  and  o\  respectively,  i.e.  wm  ~  N(0,  cr^). 
The  target  motion  is  assumed  to  be  purely  deterministic,  i.e.  there  is  no  process  noise 
involved. 


x 

Figure  1:  Trajectories  of  friendly  platforms  (Pi,  P2)  and  target  T;  z\  and  z2  represent  angular 
measurements,  ip  is  the  vergence  angle 

Assume  that  identical  angle-only  trackers  are  placed  on  each  of  the  friendly  platforms. 
Then  we  may  consider  just  one  of  them,  for  example  the  one  on  platform  1.  The  simple 
architecture  is  shown  in  Figure  2.  The  local  sensor  is  sending  angular  measurements 
z\(k)  at  every  sampling  instant  tk-  The  on-board  inertial  navigation  system  (INS)  is 
supplying  the  own-ship  state  vector1  si (k)  —  [a;i(A;),  yi{k),  yi(k)]T  at  the  same 

time.  Occasionally,  at  time  instants  ti  =  l  T0  where  l  —  b  +  T]k,  the  tracker  receives  via  the 
communication  link  additional  angular  measurements  z2(l)  and  own-ship  state  vector  s2(l) 
from  the  other  platform.  The  assumption  is  that  both  Sx (A:)  and  s 2(l)  are  known  exactly. 
Parameter  6  =  0,1,2,...  defines  the  delay  in  the  initial  measurement  transmission  while 
77  =  1,2,3,...  is  a  function  of  bandwidth.  The  case  6  =  0,  77  =  1  corresponds  to  the  full 
connection  scheme,  where  all  measurements  are  exchanged  instantaneously  between  the 
platforms.  The  case  b  >  N  represents  the  opposite  single-observer  or  no-link  situation, 
studied  for  instance  in  [4].  Cases  i]  =  2, 3,  etc.  represent  the  cases  where  only  every  second, 
third,  etc.  measurement  is  sent  over  the  link.  The  relationship  between  the  bandwidth  B 
and  parameter  ij  is  straightforward: 


VT0 

Sectors  are  represented  by  boldface  lower  case  letters. 


(2) 
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where  Nf,  is  the  total  number  of  bits  per  message2.  Given  a  link  with  a  specified  bandwidth 
B0 ,  parameter  77  can  be  calculated  as: 


77  = 


1  BnTn  1 


(3) 


Next  we  investigate  the  influence  of  the  tactical  data  link  bandwidth  on  angle-only  tracking 
accuracy. 


Figure  2:  Tracking  architecture 


3  Cramer-Rao  bound 

Taking  into  account  that  the  target  is  travelling  at  a  constant  velocity  along  a  straight 
line,  the  unknown  target  state  vector  can  be  represented  as  a  function  of  time  index 
k  =  1, . . . ,  N  as  follows: 

s(k)  =  [si(k)  s2{k)  s3(k )  54(/c)]T  =  [xt{k)  xt  yt(k)  yt]T .  (4) 

The  covariance  matrix  of  the  unbiased  estimator  s{k)  is  bounded  by  the  inverse  of  the 
Fisher  Information  matrix3  (FIM)  J  : 

E{([s(As)  -  s{k)][s(k)  -  s{k)]T}  >  J_1(fc)  (5) 

The  FIM  is  defined  as4  [9]: 

J(s)  =  E{[Vslogp(Z,s)]  [vslogp(Z,s)]T}  (6) 

where  Z  is  the  cumulative  measurement  set  which  represents  a  union 

Z  =  Z\  U  Z2.  (7) 

*A  message  from  platform  m  typically  would  consists  of  a  triple  {zm(t,);  sm(t,);  tt} 

3Matrices  are  represented  by  upper  case  boldface  letters. 

In  order  to  simplify  notation  we  drop  out  index  k  unless  there  is  possible  confusion. 
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Measurements  Z\  are  from  platform  Pi  and  measurements  Z2  are  from  platform  P2,  and 
at  time  index  k  they  can  be  expressed  as: 

Zi{k)  =  {zi(i)\i  - 

Z2{k)  -  {z2{j);j  =  1  +  b  +  rjl\  l  =  0, . . . ,  \^]  -  b  -  1}  (8) 

Density  p(Z,s)  in  (6)  is  a  joint  measurement-state  probability  density  function.  Note  that 
in  our  model  (described  in  Section  2)  both  measurements  and  target  state  are  random. 
The  target  state  s  is  random  despite  the  assumption  of  the  deterministic  evolution  of  the 
target  state  (i.e.  no  process  noise),  because  the  initial  target  state  s(0)  is  assumed  to  be 
random.  The  expectation  operator  E{  }  in  (6)  is  taken  with  respect  to  both  Z  and  s. 
Finally,  ys  in  (6)  is  the  gradient  operator  with  respect  to  s. 


3.1  The  Joint  Measurement-State  Density 

Due  to  the  independence  between  measurements5  Z\  and  Z2,  the  joint  measurement- 
state  density  function  p(Z,  s)  can  be  expressed  at  time  index  k  as  follows: 

P[Z{k),s(k)]  =  P[s{k)]-p[Z{k)\s{k)] 

=  p[s(k)]  ■p[Z1(k)\s(k)]-p[Z2(k)\S(k)}.  (9) 

Probability  density  p[s(fc)]  represents  the  prior  density  of  the  target  state  which  is  either 
known  or  can  be  worked  out  using  the  prior  information  about  the  sensor  and  the  target. 
Typically  one  knows  the  minimum  and  maximum  sensor  range,  rmjn  and  rmax,  respectively, 
as  well  as  the  maximum  target  speed  nmax.  The  common  approach  is  to  assume  a  Gaussian 
prior  density  at  k  =  0,  that  is,  s(0)  ~  A/”(so,  Po)-  The  prior  density  is  then: 

p[s(0)]  =  (2tt)^7[^[  eXP  {~^S  ~  S°^  P°  1  (s“So)}'  (10) 

The  mean  of  density  p[s(0)],  denoted  as  so,  is  either  a  priori  known  or  is  calculated  using 
the  previous  (at  k  =  0)  angular  measurement  z\  (0): 

s0  =  [f  sinzi(O)  +  xi(0),  0,  f  cos  zi(0)  +  yi(0),  0]T,  (11) 

where  f  =  (rmax  +  ^min)/2  and  si(0)  is  the  platform  1  own-ship  state  vector  at  k  =  0.  The 
covariance  of  p[s(0)]  is  typically 

P0  =  diag(of,  o2t,o2v)  (12) 

where  ar  and  ov  are  chosen  in  such  a  way  that  the  prior  density  covers  the  entire  span  of 
the  target  range  (from  rmin  to  rmax)  and  velocity  (from  -nmax  to  +nmax).  Next  we  need 
to  work  out  the  contribution  of  the  prior  density  at  time  index  k.  Note  that  due  to  the 
purely  deterministic  target  motion  the  following  linear  relationship  holds: 

s (k)  -  Ffcs(0)  (13) 

5The  independence  of  the  measurements  is  due  to  the  independence  of  the  measurement  noise  on  the 
platforms. 
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where 
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Then  it  follows  (see  for  instance  [7,  Theorem  2.11])  that  s (k)  ~  A7(sfc,  P*),  where 

Sfc  =  Ffcs0  and  P*.  =  FkP0[Fk]T. 


(14) 


(15) 


In  expression  (9),  p[Zi(k)\s(k)]  represents  the  likelihood  function  due  to  the  own-ship 
measurements.  Since  measurement  noise  w\(k)  in  (1)  is  Gaussian, 


P[zi(*)i»(*)i = n 


(TiVZjt 


exp  <  - 


[zi{i)  -  /ii(s,i)]2 


2a\ 


where  (keeping  in  mind  eq.(4)) 


hi(s,  i)  —  arctan 


sijk)  -  (fc  -  i)T0s2(k)  -  xi(i) 
s3(k)  -  (k-i)  T0  s4(k)  -  yx(i) ' 


(16) 


(17) 


Finally  the  third  density  term  on  the  right-hand  side  of  (9)  is  the  likelihood  function 
due  to  the  measurements  from  platform  2.  The  same  reasoning  as  above  leads  to  the 
following  expression: 


p[z2(k)\S(k)}=  n 


i=(l+b):ri 


exp 


\z2{i)  -  Ms,?)]2 

2cr| 


(18) 


where  notation  i  =  (l  +  b):p  signifies  that  p  is  the  increment  (i.e.  index  i  takes  the  values 
1  +  6, 1  -I-  b  +  77, 1  +  b  +  2p, . . . ).  Similarly  to  (17)  we  have 


/i2(s,  i)  =  arctan 


si(k)  -  (fc  -  i)T0S2(k)  -  x2(i) 
ss(k)  -  (k-i)  T0  s4(k)  -  y2(i) ' 


(19) 


Note.  The  expression  for  the  joint  density  function  (9)  could  be  generalised  in  a 
straightforward  manner  to  the  case  of  three  or  more  platforms  by  adding  the  appropriate 
likelihood  functions  in  the  product. 


3.2  Fisher  Information  Matrix 

Before  calculating  the  FIM  of  eq.(6)  let  us  consider  the  negative  logarithm  of  p(Z,  s), 
\(Z,s)  =  —  logp(Z,  s) 

1  ^ 

=  C  +  -(s  -  sk)T  p-1  (s  -  Sfc)  +  A-  5^[*i(t)  -  hx(s,  i )]2 

1  i=l  ,  , 

x  fc  (20) 

+  5Z  M)  -h2(s,j)}2 

2  j=l+b:r) 
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where  C  denotes  a  constant.  In  a  shortened  notation  (20)  can  be  expressed  as: 

A(Z,  s)  =  C  +  Ao(s)  +  Ai  (Zi,  s)  +  A2(Z2,  s).  (21) 

Prom  the  definition  of  the  FIM  in  (6)  we  have: 

J(s)  =  E{[Vs(C  +  A0  +  Ai  +  A2)]  [\7s{C  +  A0  +  Ai  +  A2)]T}  (22) 

Using  the  linearity  of  the  gradient  operator  and  the  mutual  independence  of  Zi,  Z2  and 
s(0),  we  can  write 

J  =  Jo  +  Ji  +  J2  (23) 

where 

Jm  —  E{[Vs^m]  [Vs^m]  }  (m  =  0, 1,  2) 

Next  we  derive  each  term  in  (23). 


3.2.1  Derivation  of  J0. 


Note  that  [6] 


„  X  M  _  ldKs-SfcfPfc^s-s*)] 

VsAo(s)  -  - - ^ - 


* k  (s  -  sk )  (24) 


Then  using  the  basic  rules  of  matrix  algebra  and  the  fact  that  the  inverse  of  the  covariance 
matrix  is  a  symmetric  matrix  we  have 


J„  =  E{Pp(s-sl)((s-s*)T[Pjr1]7'} 

=  P*~‘E{(s  -  s*)((s  -  stflPj 1 

=  (25) 


3.2.2  Derivation  of  Ji  and  J2. 


From  [4]  it  follows  that: 


and 


where 


Ji 


=  E 


^[VsMM]  [v»MM)]r  f 

i=  1  J 


J2  =  E 


1  ' 

-2  5Z  [VsMM)j  [Vs^2(s,x)]T 

a2  i=l+6:r) 


Vs^m 


dhm 

dsi 


dhm  dhm 
ds2  dss 


{m  =  1,2), 


(26) 


(27) 


(28) 


6 


DSTO-TR-1053 


and  with  expectation  E  in  eqs.(26)  and  (27)  taken  over  the  prior  distribution  of  s.  The 
appropriate  terms  in  (28)  are  as  follows: 


dhm 

dsi 

dhm 

ds2 

dhm 

ds3 

dhm 

ds4 


yt(i)  -  ym(i) 


[*t(0  -  Xm{i)]2  +  [ yt{i )  -  ym(i)}2 

,  ,,  dh  777 

—  {k  -  i)T0 


ds  i 

xt{i)  -  xm{i) 


\xt{i)  -  xm(i)}2  +  [ yt(i )  -  ym(i )]2 
dhm 


(k  -  i)T0 


ds3 


Let  us  introduce  the  following  notation  for  the  relative  state  vectors: 


—  s  8777  —  [xm,  xm,  ym,  ym] 


( m  =  1,2) 


(29) 


Note  that  sm(A;)  ~  jV(sfc  —  sm(fc),Pfc)  where  s^.  and  Pfc  were  given  by  (15).  It  is  now 
straightforward  to  verify  that  in  order  to  evaluate  the  4x4  matrices  Ji  and  J2  we  need 
to  work  out  the  expressions  for 

f  1  r  xy 

E  l  (52  +  y2)2  J  {mdE\W+W 

where  x  and  y  are  both  Gaussian  and  independent.  This  step  can  be  performed  by 
numerical  integration  or  Monte  Carlo  averaging  [8]. 

The  FIM  given  by  (23)  depends  on  (i)  the  geometry  of  the  scenario  (the  motion  of  the 
communicating  platforms  and  the  target);  (ii)  measurements  accuracy  (<xi  and  a2)]  (iii) 
sampling  interval  (T0);  (iv)  communication  bandwidth  (77).  The  CRLBs  are  calculated  as 
the  diagonal  elements  of  the  inverse  FIM,  i.e. 


CRLB{%}  -  [J-1]jj  ( j  =  1,2, 3, 4)  (30) 

where  Sj,  j  =  1,2, 3, 4  are  expressed  in  the  Cartesian  coordinates,  see  eq.(4).  In  practice, 
however,  it  is  often  of  more  interest  to  consider  the  estimate  of  target  range  and  its  CR 
lower  bound.  This  CRLB  will  be  derived  next. 


3.3  CRLB  for  the  Target  Range  Estimate 

The  target  range  is  a  non-linear  function  of  the  state  vector  s: 

Rt(k)  =  \/x2t(k)  +  y?{k)  =g[s{k)}  (31) 

According  to  [9,  pg-83] 

<«> 

i  j  J 
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The  partial  derivatives  are  as  follows: 


dg 

dsi 

dg 

ds2 

dg 

dsz 

dg 

ds4 


which  leads  to  the  following  expression: 


(33) 

(34) 

(35) 

(36) 


CRLB(Rt) 


_|_CRLB(it)  +  _^_ 


+  (J-'U  + 


?t^crlb®,) 

xt  '  Vt 

(37) 


Note.  If  estimates  x t  and  yt  are  efficient  (meet  their  respective  CRLBs)  the  estimate 
of  range  Rt  =  y/xf  +  y\  may  not  be  efficient  because  g(s)  is  non-linear  [9,  pg.84j. 

The  next  section  will  discuss  the  influence  of  bandwidth  on  the  tracking  performance 
and  in  addition  will  compare  two  algorithms  for  two-platform  angle  only-tracking  with  the 
derived  CRLB. 


4  The  Influence  of  Bandwidth 

4.1  Discussion 

Consider  the  scenario  shown  in  Figure  3  for  the  discussion  that  follows.  A  pair  of 
communicating  platforms  flies  in  parallel  at  a  speed  of  540  km/h  with  a  mutual  distance 
of  approximately  6  km  during  the  first  half  of  the  scenario  (k  =  1, . . . ,  30).  Then  platform 
1  performs  a  sharp  turn  in  order  to  increase  the  vergence  angle,  and  thereafter  both 
platforms  again  move  at  the  constant  speed  and  heading.  Target  T  is  moving  eastwards 
at  a  speed  of  790  km/h.  The  sampling  interval  is  T0  =  2  sec. 

Since  the  considered  scenario  (including  the  initial  target  state  s(0))  is  purely  deter¬ 
ministic,  the  calculation  of  the  CRLBs  can  be  somewhat  simplified.  The  uncertainty  about 
s(0)  can  still  be  expressed  in  terms  of  the  prior  covariance  matrix  Po,  but  there  is  no  need 
to  perform  the  expectation  operation  in  eqs.(26)  and  (27).  For  the  purpose  of  studying 
the  influence  of  bandwidth  this  will  be  adequate. 

The  CRLBs  of  the  range  estimates  as  a  function  of  time  index  k  are  shown  in  Figure  4 
for  b  =  0  and  r\  =  1, 2,4,8, 16  and  no-link  (single  Pi  observer)  case.  The  prior  covariance 
matrix  Po  of  (12)  is  chosen  with  parameters 

ar  =  25  km  and  av  =  1080  km/h.  (38) 

The  assumed  standard  deviation  of  angular  measurements  is  o\  =  a2  =  0.5°.  Consider 
the  first  extreme,  the  ‘no-link’  case.  Note  that  at  k  —  1  the  CRLB  is  determined  by  the 
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Figure  3:  Angle-Only  Tracking  scenario 


prior  covariance  matrix  (the  square  root  of  the  CRLB  corresponds  roughly  to  ar).  As 
the  time  progresses  the  CRLB  varies  slightly  due  to  the  geometry  and  the  collected  angu¬ 
lar  measurements  until  the  manoeuvre,  when  it  sharply  decreases  as  the  target  becomes 
“observable”.  The  opposite  extreme  is  the  ‘full  connection’  case  (q  =  1,  6  =  0).  In  this 
case  the  CRLB  is  significantly  lower  than  in  the  ‘no-link’  case  and  completely  indepen¬ 
dent  of  the  manoeuvre  as  the  target  is  observable  throughout  the  scenario.  The  cases 
q  —  2,4, 8, 16  are  between  the  two  extremes,  and  the  respective  CRLBs  follow  the  pattern 
of  the  received  measurements  over  the  communication  link. 

Consider  for  example  the  CRLB  curve  for  q  =  16  and  6  =  0  case  in  Figure  4.  Af¬ 
ter  receiving  the  first  angle  measurement  from  the  other  platform  (at  k  =  1)  the  range 
estimate  error  standard  deviation  is  about  12  km.  Due  to  the  lack  of  measurements  for 
16  consecutive  time  steps,  the  standard  deviation  will  rise  to  the  level  which  corresponds 
to  the  ‘no- link’  case  (about  25  km).  At  k  =  17  when  the  link  becomes  active  again,  the 
range  error  standard  deviation  drops  dramatically  ,  although  not  quite  to  the  level  which 
corresponds  to  the  ‘full- connection’  case.  Then  it  continues  to  rise  again  until  the  link 
is  active  again,  etc.  Interestingly,  after  the  manoeuvre  practically  all  q  =  2, 4, 8, 16  cases 
converge  to  the  full  connection  performance  (after  k  =  50  their  CRLBs  are  similar).  For 
6  >  0  the  CRLBs  are  identical  to  the  ‘no- link’  case  for  1  <  k  <  1  +  6.  During  the  remaining 
interval  l  +  b  <  k  <  N,  the  CRLBs  are  identical  to  the  6  =  0  cases  shown  in  Figure  4. 

Figure  5  repeats  the  CRLBs  of  the  range  estimate  errors,  but  this  time  the  assumed 
measurement  accuracy  is  5  times  better,  i.e.  <7i  =  02  =  0.1°.  Notice  that  with  the  higher 
quality  of  angular  measurements  the  need  for  a  high  capacity  link  is  somewhat  reduced, 
and  the  cases  q  =  2  or  even  q  =  4  are  almost  as  good  as  the  full  connection  case. 


4.2  Tracking  Algorithms  Against  the  CRLB 


For  completeness  this  section  investigates  the  performance  of  a  pair  of  algorithms  for 
two-platform  angle  only  tracking  and  compares  their  performance  with  the  CRLBs  derived 
and  discussed  earlier. 
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xIO4 


Figure  4:  Cramer-Rao  lower  bounds  of  range  estimate  standard  deviation  for  01  =  cr?  =  0.5°. 
The  link  parameters:  6  =  0  and  t)  =  1,2, 4,  8, 16  and  no-link  case. 

4.2.1  MLE  Algorithm 

The  first  algorithm  we  consider  is  the  Maximum  Likelihood  estimation  (MLE)  algo¬ 
rithm.  We  use  MLE  to  estimate  the  range,  keeping  in  mind  the  note  made  at  the  end  of 
Section  3.3.  The  MLE  is  a  batch  algorithm  although  we  are  interested  in  the  performance 
as  a  function  of  time  index  k.  Hence  we  apply  the  MLE  progressively  at  each  time  step, 
over  a  cumulative  measurement  set.  The  MLE  implemented  in  this  way  is  very  close  to 
an  efficient  estimator  of  range.  Note  however  that  MLE  does  not  use  prior  information 
and  hence  is  appropriate  only  for  the  “observable  target”  segments  and  larger  values  of  k 
where  the  significance  of  prior  information  diminishes  anyway.  In  mathematical  terms  the 
MLE  is  expressed  as 


sMle  =  argminA(s)  (39) 

S 

where  A(s)  is  given  by  (20).  Minimisation  was  performed  as  in  [4]  using  a  MATLAB© 
routine  fmins. 


4.2.2  EKF  with  Triangulations 
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Next  we  consider  a  recursive  algorithm  which  uses  the  prior  information:  an  Extended 
Kalman  filter  (EKF)  in  the  Cartesian  coordinates  with  the  use  of  triangulated  range 
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x  104 


Figure  5:  Cramer-Rao  lower  bounds  of  range  estimate  standard  deviation  for  o\  —  02  =  0.1°. 
The  link  parameters:  6  =  0  and  77  =  1, 2, 4, 8, 16  and  no-link  case. 


estimates.  This  EKF  is  by  the  assumed  convention  (in  relation  to  Figure  2.)  installed  on 
platform  1,  and  its  measurement  vector  is  defined  as  follows: 


z(k)  = 


z\ 

{Ri,  zi)T 


iik^l, 
iik  =  l. 


(40) 


Equation  (40)  indicates  that  when  angular  measurements  from  two  platforms  are  available 
(if  k  =  l)  the  measurement  vector  is  augmented  by  the  target  range  estimate  (from  platform 
1)  obtained  by  triangulation  [3,  Sec.  10.6.4], 


Ri 


d  ■  sin  62 
sin 'ip 


(41) 


In  order  to  explain  eq.(41)  we  introduce  unit  vectors  ui  and  112  defined  as  the  directions 
from  platform  Pi  and  P2  respectively  to  target  T  (see  Figure  1).  They  are  calculated  from 
measurements  z\  and  22  as  follows: 


um  =  sin  2^  +  cos  zmiy  (m  =  1, 2) 

where  ix  and  iy  are  unit  vectors  in  the  directions  of  the  x  and  y  axes,  respectively.  We 
also  need  to  introduce  the  baseline  vector  d  defined  by  the  own-ship  states  Si  and  S2 
as:  d  =  (x2  —  Xi)ix  +  (j/2  —  yi)iy-  Now  the  quantities  in  (41)  are  as  follows:  d  =  |d|, 
ip  =  arcsin(|ui  x  U2I)  and  02  =  arcsin(|d  x  U2I/0J).  The  variance  of  the  triangulated  range 
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is  given  by  [3,  Sec.  10.6.4] 


R\  cos2  ip  aji  +  /?2<7 


6  2 


R  i 


sin  ip 


(42) 


where  R2  =  d  ■  sin#i/ sin  ip,  is  the  triangulated  range  from  platform  P2.  Since  own-ship 
state  vectors  are  known  exactly,  a6y  =  o\  and  ag2  =  cr2.  Introducing  the  relative  target 
state  vector  (with  respect  to  platform  1)  as  in  eq.(29), 

s{k)  =  s (k)  -  si (k)  =  [Si,  h,  «3i  «4]T 


where  s(fc)  is  the  target  state  vector  of  eq.(4),  the  state  equation  can  be  expressed  as 
follows: 


s  (k  +  1)  =  F  s  (k)  -  U(fc  +  1,  k) 
where  F  was  defined  by  (14)  and 


(43) 


u(fc  +  i,*o 


xi (k  +  1)  -  Xi (k)  -  T0&i(k ) 
±i(k  -f  1)  —  i\ (k) 
yi(k  +  l)  -yi(k)  -Tayi{k) 
yi{k  +  1)  -  yi{k) 


is  a  vector  of  deterministic  inputs  which  account  for  the  observer  (platform  1)  acceleration. 
The  measurement  equation  is  non-linear, 


z  (k)  —  fi[s(fc)]  +  w  (k). 

The  non-linear  function  h[  ]  in  (44)  is  defined  as6: 


h[  s]  = 


VH  +  S3 

arctan  41 


(44) 


(45) 


The  covariance  matrix  of  w (k)  in  (44)  is  R  =  diag[a^,  a2],  where  o2^  is  given  by  (42). 
The  standard  EKF  equations  [2,  Sec. 10.3.3]  are  applied  with  the  equivalent  measurement 
matrix  given  by: 


Hfc+x  = 


dh[s(k  +  1)] 
ds 


s=s(A-+ l|fc) 


a2  -2 
S1+S3 


--2  ,  :2 
*l+«3 


0 

0 


0 

0 


£+llfc 


(46) 


4.2.3  Results 

The  simulations  were  performed  for  the  scenario  described  in  Section  4.1.  The  simu¬ 
lation  results  are  shown  in  Figures  6. (a),  (b)  and  (c)  for  cases  (i)  no-link  (ii)  r]  =  8,  b  =  0 
and  (iii)  ij  =  1,  b  =  0  (full  connection),  respectively.  In  all  cases  the  measurement  noise 
was  set  to  a\  =  a2  =  0.1°,  and  the  sampling  interval  was  T0  =  2  sec.  The  y-axis  in  all 

6 We  consider  only  k  =  l  case  of  (40),  since  the  other  one  follows  directly. 
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the  figures  represents  the  root-mean-square  (RMS)  range  errors.  The  calculated  CRLBs 
are  always  shown  with  dashed  lines.  The  MLE  curves  are  obtained  from  simulations  by 
averaging  over  100  independent  Monte  Carlo  runs.  As  expected  they  match  the  calculated 
CRLBs.  The  EKF  algorithm  uses  the  prior  knowledge  about  the  target  state  in  the  form 
of  the  prior  density  (10).  In  the  simulations  the  mean  of  the  prior  density  was  set  to  the 
true  value  of  s(0),  while  the  covariance  was  the  one  used  in  the  calculations  of  CRLBs 
[see  (38)].  The  EKF  curves  are  obtained  by  averaging  over  500  Monte  Carlo  runs,  and  are 
shown  with  dash-dotted  lines  with  triangles.  Since  the  EKF  is  a  suboptimal  algorithm,  its 
RMS  errors  are  always  higher  than  the  calculated  CRLBs,  although  note  that  in  all  the 
cases  considered  they  follow  the  trend  of  the  respective  bounds. 


5  Summary 

The  influence  of  communication  link  bandwidth  on  angle-only  tracking  with  a  pair 
of  communicating  moving  platforms  is  studied.  The  method  of  investigation  is  based  on 
the  Cramer-Rao  lower  bounds  which  are  derived  for  recursive  estimators  with  prior  in¬ 
formation.  The  bounds  were  confirmed  by  Monte  Carlo  simulations  using  the  maximum 
likelihood  estimator  over  a  cumulative  measurement  set  and  the  extended  Kalman  filter 
with  triangulation.  The  tracking  performance  in  the  limited  bandwidth  environment  is 
shown  to  be  between  the  two  extremes,  the  full-connection  and  the  no-connection  cases, 
and  it  follows  the  pattern  of  the  received  measurements  over  the  link.  As  the  time  pro¬ 
gresses  the  limited  (but  non-zero)  bandwidth  cases  tend  to  approach  the  full  connection 
performance.  This  effect  is  more  pronounced  when  angular  accuracy  is  higher. 

The  study  presented  in  this  report  was  carried  out  under  a  number  of  assumptions 
about:  the  target  (non- manoeuvring  single  target),  the  sensor  (no  false  detections  and  no 
missed  detections)  and  the  data  link  (no  transmission  delay).  The  future  work  will  try  to 
relax  some  of  these  assumptions. 


13 


s 


B 


Figure  6:  The  RMS  range  errors  of  the  MLE  and  the  EKF+triangu 
Simulations  results  for:  (a)  no-link  case;  (b)  r\  =  8,  b  =  0  case 
connection)  case. 
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